Thallium 7p lifetimes derived from experimental data 
and ab-initio calculations of scalar polarizabilities 
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Two different theoretical methods have been used to complete a new calculation of polarizability 
in the thallium 6pi/2, 7s, and 7pi/2 states. The predictions of the two methods agree to within 1% 
for the 6pi/2 and 7s states and 2% for 7p]/ 2 state. We find that the theoretical expression for the 
6P1/2 ^7s transition polarizability difference, Aao, is dominated (greater than 90% contribution) by 
mixing of the 7s state with the 7p!/ 2 and 7p 3 / 2 states. By comparing the theoretical expression to an 
existing measurement of Aao [Doret et ai, Phys. Rev. A 66, 052504 (2002)], new, highly-accurate 
values for the thallium 7p excited-state lifetimes have been extracted. The scalar polarizability of 
the 7pi/2 state is also computed, anticipating an experimental determination of this quantity, which 
will then enable a high-precision determination of the 6dj — 7pj< transition rates and provide a 
benchmark test of the two theoretical approaches in the near future. 



I. INTRODUCTION 

Thallium has played an important role in atomic- 
physics-based tests of discrete symmetry violation over 
recent decades [IHl]. The size of these symmetry- 
violating observables scales rapidly with the atomic num- 
ber, encouraging the use of high-Z systems. This there- 
fore requires independent, precise atomic wave function 
calculations in order to distinguish quantum mechanical 
effects from the elementary particle physics observables 
being targeted. For example, theoretical uncertainties in 
ah initio wave function calculations in thallium currently 
limit the quality of the standard model test provided by a 
1995 thallium parity nonconservation measurement 0, 0| ■ 
Similarly, high-precision atomic theory is essential to 
interpret results from searches for atomic and molecu- 
lar electric dipole moments (EDMs), as evidenced by 
recent calculations of the thallium EDM enhancement 
factor @, [B-0|- Independent, high-precision atomic struc- 
ture measurements serve as an important tool in testing 
the accuracy and guiding the refinement of theoretical 
techniques for multi-electron systems such as thallium. 
Over the recent years, we have completed precise mea- 
surements of thallium transition amplitudes [8], hyperfinc 
splittings [9], and polarizab ility (Toj which show excellent 
agreement with theory [TTJ, Il2j . More recently, a similar 
theoretical approach to that used for thallium has also 
been applied to other trivalent Group IIIA systems such 
as indium and gallium [ill [l4j|. 

Very recently, theoretical and experimental work have 
come together in the indium atomic system. A mea- 
surement of the Stark shift within the indium 410 nm 
5pi/2 ~ 6S1/2 transition l_5j yielded a value for the 6s — 
5pi/2 polarizability difference with 0.3% uncertainty. At 
the same time, a new ah initio theory effort, using two 
complementary, high-precision techniques, yielded a the- 
oretical value for this quantity in excellent agreement 
with the experimental result, and with 2% estimated 



uncertainty[3|. Because the theoretical expression for 
the 6s polarizability is dominated by terms involving the 
6s — 6pi/2 and 6s — 6^3/2 mixing, we show in Refs. [lilllsT] 
that a comparison of experimental and theoretical results 
can produce new values for the 6p-state lifetimes with un- 
certainties below 1%. 

In this paper, a similar approach of combining high- 
precision calculation and experiment is applied to thal- 
lium. An extensive calculation including uncertainties 
is undertaken using both a coupled-cluster (CC) as well 
as a configuration interaction + all-order (CI+all) ap- 
proach to compute the polarizability of the thallium Qp\/2 
ground state, as well as the 7s and 7p\/2 excited states. 
We use a comparison of this theory to the 2002 thallium 
Stark shift result [l(| to extract the most precise values to 
date for the thallium 7p-state lifetimes. We also outline 
ongoing experimental work which will allow precise mea- 
surements of excited-state Stark shifts in both thallium 
and indium. Such measurements will then be combined 
with theoretical polarizability results to accurate predict 
the thallium 6d— 7p and indium 5d — 6p transition rates. 



II. CALCULATION OF POLARIZABILITIES 

The valence static polarizabilities of a Tl atom can be 
calculated as sum over states: 
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where the sum over n runs over all states with allowed 
(J||D||J n ) electric-dipole transitions. 

The 7s — 6pi/2 Stark shift in Tl (i.e. the difference of 
the 7s and &P1/2 polarizabilities) is strongly dominated 
by the contributions from the 7s — 7pj transitions to the 



7s polarizability. Therefore, 
this Stark shift carried out in 



accurate measurement of 
[loj can be used to extract 
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TABLE I: Contributions to the 7s, 6P1/2, an d 7pi/2 static 
polarizabilities are given in Oq in columns labeled "«o" . 
The experimental energies [161 ] (in cm -1 ) and the theoreti- 
cal electric-dipole reduced matrix elements (in a.u.) used to 
calculate dominant contributions are listed in columns labeled 
"A_B" and "D" . The CC and Cl+all-order electric-dipole ma- 
trix elements and the polarizability contributions are listed in 
columns labeled "CC" and "CI-I-A11", respectively. 



Contribution 


AE 
Expt. 


CC 


D 
CI+A11 


CC 


CI+A11 


7s polarizability 


6f>l/2 
7f>l/2 

8P1/2 

(n > 8)p 1/2 


-26478 
7682 
14891 


1.826 
6.016 
0.706 


1.798 
6.050 
0.693 


-9.2(4) 
345(12) 
2.4(5) 
0.7(2) 


-8.9 
349 
2.4 


6P3/2 

7p 3 /2 

8P3/2 

(n > 8)p 1/2 


-18685 
8684 
15263 


3.397 
8.063 
1.474 


3.395 
8.108 
1.509 


-45(2) 
548(22) 
10(1) 
5(1) 


-45 
554 
11 


Core 
Total 










24(1) 
881(25) 


5 

887 






6p 


x/2 polarizability 






7s 

8s 

(n > 8 


)* 


26478 
38746 


1.826 
0.535 


1.798 


9.2(5) 
0.54(5) 
0.8(3) 


8.9 


6^3/2 

7d 3/2 
(n > 8 


) rf 3/2 


36118 
42011 


2.334 
1.101 


2.377 


11.0(4) 
2.1(1) 
6.4(2.8) 


11.5 


Core 
vc 

Total 










24.1(1.2) 
-4.2(9) 
50.0(3.0) 


5.0 
-0.4 
50.7 






7 P 


1/2 polarizability 






7s 

8s 

(n >8 


)s 


-7682 
4586 


6.013 
6.189 


6.050 


-344(3) 
611(5) 
22(1) 


-349 


6^3/2 

7d 3 /2 

(n > 8)d 3/2 


1958 
7852 


10.726 
4.767 


10.649 


4298(24) 
212(13) 
95(13) 


4237 


Core 
Total 
Final 










24(1) 
4918(30) 
4918(120) 


5 

4831 



7 s — 7pj matrix elements if all other smaller contributions 
to the 6pi/2 and 7s polarizabilities are calculated. The 
extraction of the matrix elements also require the eval- 
uation of relevant theoretical uncertainties. Combining 
7s — 7pj matrix elements with experimental transition 
energies gives the 7pi/ 2 and 7p 3 / 2 lifetimes. 

The polarizabilities of the ground 6pi/2 and excited 7s 
states were calculated using the linearized coupled- clus- 
ter method in [121 ] , but their uncertainties were not eval- 
uated. In this work, we carry out additional linearized 
coupled-cluster calculations (CC) to estimate the uncer- 
tainties of each term. We also carry out another inde- 
pendent calculation of these polarizabilities using a re- 
cently developed hybrid approach that combines configu- 
ration interaction and linearized coupled-cluster method 



(Cl+all-order) [T{|. This calculation treats Tl as a sys- 
tem with three valence electrons and accurately accounts 
for configuration mixing and valence- valence correlations. 
The CC calculation treats Tl monovalent system 

with 6s 2 considered to be a part of the closed core; how- 
ever, the CC approach includes some additional high- 
order corrections to the dipole operator. Since these two 
methods differ in their inclusion of higher-order effects, 
comparing their results provides additional evaluation of 
the uncertainty of our calculations. We refer the reader 
to Refs. [HMmi and HGIHHH for detailed descrip- 
tions of the linearized coupled-cluster and Cl+all-order 
methods, respectively. The results of both methods were 
recently compared for In polarizabilities in [l4j . 

The breakdown of the contributions to the 7s and 6pi/2 
Tl polarizabilities is listed in Table [I] Experimental en- 
ergies are given for all terms that are listed separately, 
such as 7p contributions to the 7s polarizability. The 
uncertainties of the main CC terms are determined from 
the spread of four different coupled-cluster calculations 
carried out in this work (with and without the pertur- 
bative triple terms and with inclusion of the scaling to 
account for some missing higher-order corrections); CC 
matrix elements from [l2j are kept as final values. The 
determination of the uncertainties is described in detail 
in Refs. [13, EJ. The (n > 8) contributions for the 7s po- 
larizability contain 9p contribution calculated using the 
all-order method and all other (n > 9) terms calculated 
in the random-phase approximation (RPA) and scaled to 
account for higher-order corrections. The scaling factor is 
determined as the ratio of the total CC value for the main 
n = 6 — 9 terms and corresponding RPA result. For the 
6pi/2 polarizability, (7— 10)s and (6— 9)d 3 / 2 contributions 
are calculated by combining CC matrix elements and ex- 
perimental energies, and the remaining contributions are 
calculated together using the scaled RPA approach. The 
difference of the ab initio RPA and scaled RPA values is 
taken to be uncertainty of these high— n contributions. 
The ionic core polarizabilities and small (vc) term that 
accounts for the occupied valence shell(s) are listed sepa- 
rately in rows "core" and "vc" . The vc term is negligible 
for the 7s and 7px/2 polarizabilities. Core and vc contri- 
butions are calculated in the RPA. The differences of the 
Dirac-Fock and RPA values are taken to be their uncer- 
tainties. We note that core polarizability is much larger 
in the CC method, since 6s 2 shell is included in the core 
in the CC calculation, while the 6s shell belongs to the 
valence space in the trivalent Cl+all-order calculation. 

The sum-over-states is not used in the Cl+all-order 
calculation of the polarizabilities which is carried out by 
solving the inhomogeneous equation of perturbation the- 
ory in the valence space [261 ] . However, we evaluated 
a few dominant terms separately by combining Cl+all- 
order matrix elements with experimental energies to com- 
pare these terms in both approaches. These results are 
listed in the last column of Table |H 

We also calculated the thallium 7pi/ 2 polarizability 
using both CC and Cl+all-order methods. Since the 
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TABLE II: Final values of the 7s and Qpi/2 polarizabilities and their difference Aao (a.u.). Determination of the reduced 
electric-dipole 7s — 7pj matrix elements (in a.u.) and 7pj lifetimes (in ns) from the combination of measured Stark shift [13] 
and theoretical values. The quantity C is the value of Aao(7s — 6^1/2) with the contribution of the 7s — 7pj transitions 
subtracted out. a Ref. [HI, 6 Ref. [hJ, c Ref. [|, d Ref. @], e Ref. [l]|. 





Qo(7s) 


oo(6pi/ 2 ) 


Aq (7s - 6pi /2 ) 


C 


D(7s - 7 Pl/2 ) 


D(7s - 7p 3 / 2 ) 


r(7pi/ 2 ) 


t(7p 3/2 ) 


cc 


881 


50.0 


831 


-61.7 


6.016 


8.063 


60.15 


46.38 


CI+AU 


887 


50.7 


836 


-66.4 


6.050 


8.108 


59.48 


45.86 


Final 


881(9) 


50.0(1.0) 


831(8) 


-61.7(6.7) 


6.013(27) 


8.058(37) 


60.21(55) 


46.44(42) 


Expt. 

Theory 

Theory 




51(7)" 
49.2 c ,48.8 d 
52.1(1.6) e 


829.7(3.1) t 













7pi/2 — 6c?3/2 matrix element strongly dominates 7pi/2 
polarizability, our calculation can be used to extract this 
matrix element if either the 7pi/ 2 — Qpi/2 or TP1/2 — 7 s 
Stark shift is measured with high-precision. We have 
used the 7s—7pi/2 matrix element determined in the next 
section to provide more accurate recommended value. 
We determine the contribution of all other terms except 
the 7pi/2 — 6d 3 / 2 term to be 620(36) a.u. (see Table IJ). 
The determination of the final uncertainties is described 
in the next section. 



III. DETERMINATION OF THALLIUM 7p 
LIFETIMES 
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Separating the 7s — 7pj contributions (see Eq.([T])), we 
write the Aao(7s — Stark shift as 



Aa {7s-6 Pl/2 )=BS + a 



where 



B 



3 \E(7p y2 ) - E(7s) E(7 P3/2 ) - E(7s) 



(2) 



(3) 



S = D 2 is the 7s—7p 1 / 2 line strength, R is the ratio of the 
D(7s — 7p 3 / 2 ) and D(7s — 7p\/ 2 ) reduced El matrix ele- 
ments R = 1.340(4), and term C contains all other con- 
tributions to the Stark shift. We use the Aao (7s — GP1/2) 
measured in Ref. [Ioj |. Combining experimental energies 
from [l(| and our theoretical value of the ratio gives 
B = 24.65(9) a.u. The final results of our calculations 
for 6pi/ 2 , 7s polarizabilities, their difference, and term 
C are given in Table |TTJ Our theoretical value for the 
7s — 6pi/2 Stark shift is in excellent agreement with the 
experiment [lol ] . The ground state polarizability is com- 
pared with other theory 0, [|| [l8[ and experiment [ljj ■ 
Since CC and CI+all-order values include all dominant 
correlation corrections between these two calculations, 
we estimate the uncertainty in the dominant contribu- 
tions as the difference between the CC and CI+all-order 
values 5a. Then, we assume that all other uncertainties 
do not exceed the uncertainty of the dominant correc- 
tions 5a. Adding these uncertainties in quadrature, we 
arrive at he total uncertainty of \[25a. The uncertainty 



FIG. 1: (Color online) A comparison of the present 7p-state 
lifetime determinations to older experimental results. 7pi/2 
results are shown in the top portion of the figure, and 7p3/ 2 
results below. References as follows: (red) triangles: Ref. 
[27] ]; (blue) squares: Ref. [3i (g reen ) circles: present work. 



in the term C is determined by the same procedure. The 
uncertainty in the 7s—7p\/ 2 line strength S is determined 

as 



(<5Aa ) 2 + {S5Bf 



(4) 



The lifetimes of the 7pj states are obtained using r a 
1/A a b , where the transition rate A a b is given by 



A, 



2.02613 x 10 18 S ab 



A ab 



2J a + l 



(5) 



the transition wavelength X ao is in A . The final recom- 
mended values for the 7pi/ 2 and 7p 3 / 2 lifetimes are listed 
in Table |TTJ The purely theoretical lifetime results of the 
CC calculation are in very good agreement with these rec- 
ommended values. Fig. Q] compares the present results 
to earlier experimental determinations of these lifetimes. 
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IV. MEASUREMENT OF EXCITED-STATE 
STARK SHIFTS IN THALLIUM AND INDIUM 



As described here, as well as in [14j . calculations of 
the excited p-state polarizabilities, coupled with future 
measurements of excited-state Stark shifts, will allow im- 
portant new tests of the atomic theory. In particular, a 
measurement of the Stark shifts involving 7p states in 
Tl and 6p state in In will allow a definitive test of the 
CC and CI+all-order theoretical methods. While both 
approaches give results in very close agreement for the 
ground state and the first excited ns state polarizabilities 
of In and Tl, the differences increase for the next excited 
np state polarizability (to 2% for Tl and 4% for In [lij]). 
Precise Stark shift measurements involving these excited 
np states will directly address the question of whether ac- 
curate treatment of the configuration mixing or higher- 
order corrections to the matrix elements are more im- 
portant for such states. Moreover, measurement of these 
Stark shifts will allow determination of the 7p — 6d tran- 
sition rates in Tl and 6p — 5d transition rates in In. These 
are very important for improved theoretical descriptions 
of the d-state properties. These rates cannot be deter- 
mined accurately from nei-state lifetime measurements 
such as reported in [29j owing to very small branching 
ratios. 

We are currently undertaking atomic beam-based ex- 
periments to measure the excited p-state polarizabilities 
in both thallium and indium. To achieve this, having 
completed Stark shift experiments in the ground state 
410 nm and 378 nm transitions of these two Group 
IIIA elements, we are now introducing a two-step, two- 
color spectroscopy measurement scheme for use with our 
atomic beam apparatus. For both atomic systems, we 
have completed such two-step spectroscopy experiments 
in a vapor cell environment to study excited-state hyper- 
fine structure and isotope shifts [30, [H|. I n these two- 
step experiments, we begin by locking the blue or UV 
laser to the first transition step, using a supplemen tary 
atomic vapor cell and technique developed recently [32)]. 
We then overlap this laser beam with that of a second 
infrared laser and intersect both with our atomic beam 
using a transverse geometry. Interestingly, the thallium 
7s— 7pi/2, and the indium 6s— 6px/2 and 6s— 6^3/2 transi- 
tions, with resonance wavelengths of 1301 nm, 1343 nm, 
and 1291 nm respectively, can all be reached by a sin- 
gle external cavity diode laser which is currently in use 
(Sacher Lasertechnik, model TEC-150-1300-050). Using 
the FM spectroscopy technique described in pj| , will will 
extract high resolution spectra from the atomic beam 
transmission signal of the infrared laser, expected to be 
an order of magnitude weaker than the signal from the 
analogous single-step experiment, given that here we pro- 
mote only a fraction of the ground state atoms to the 
intermediate state. 

In general, extracting Stark shift information in a two- 
step transition experiment is complicated by the fact that 
both transitions are shifted when the static field is turned 



(n+1)s 




(n+1)p, 



FIG. 2: A schematic diagram showing the configuration of 
low-lying energy levels of both thallium and indium. The 
dotted lines suggest the Stark-shifted levels, where the size of 
the level shifts (not to scale) become successively larger for 
the higher lying states. 



on. Fig. [2] indicates the general form of such a three- 
level, two-step process in either indium or thallium. As 
indicated, in the presence of an electric field, all energy 
levels are Stark-shifted downwards with the magnitude 
of the shift increasing for higher-lying states. We define 
the magnitude of the Stark shift within the lower (upper) 
transition as A21 (A32), and the field-free resonance fre- 
quencies of these transitions as f 21 and f® 2 respectively. 
We note that for both elements A 32 3> A 2 i. By keeping 
the first-step laser locked to the atomic transition in a 
field-free region, the first-step excitation in the electric 
field region would then be shifted slightly out of reso- 
nance. This results in the excitation of a non-zero veloc- 
ity class of atoms. In this case, the resonance frequency 
for the second step transition, /g 2 , in the presence of the 
electric field is given by: 



/32 - /32 - 1^32 1 



|A 2 i|^, 
J21 



(G) 



where the final term results from the Doppler shift pro- 
duced by the off-resonant first-step excitation. For elec- 
tric field of 10 kV/cm, typical for these experiments, 
A 2 i « 10 MHz, still much less than the ~100 MHz resid- 
ual Doppler width in the atomic beam geometry. Thus, 
the decrease in excitation efficiency for the first step tran- 
sition in the presence of the electric field will not be signif- 
icant. Furthermore, since the optical resonance frequen- 
cies /21 and /32 in Eq. [6] are known, and A21 has been 
previously measured to high accuracy (lOj, we will be able 
to determine A32 unambiguously. For transitions involv- 
ing J = 1/2 states, there is only a scalar component to 
the polarizability, but for the indium 1291 nm 6s — 6^3/2 
transition, there will exist both scalar and tensor contri- 
butions, both of which were computed in pj]]. In this 
case we will study each resolved hyperfine transition and 
will vary the laser polarization relative to the static field 
direction to extract both polarizability components. 
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V. CONCLUSIONS 

In conclusion, through a comparison of an existing 
thallium Stark shift measurement and new ab initio cal- 
culations of scalar polarizabilities, we have derived new, 
highly- accurate values for the thallium 7p excited-state 
lifetimes. Recent measurements of Stark shifts in thal- 
lium and indium serve an important function as bench- 
mark tests of two distinct atomic theory techniques that 
can be applied to these multivalence systems. Future 
measurements of excited-state Stark shifts in both el- 
ements will test the atomic theory approaches in im- 
portant new ways since the relevant polarizability is 
dominated by mixing with excited d-states in these sys- 
tems, whose theoretical contributions are more uncertain. 



A theory-experiment comparison will thus allow precise 
derivation of thallium 7p ~ 6d and indium 6p — 5d tran- 
sition matrix elements. 
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